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Abstract We investigate tiie variation of bar strength 
with central velocity dispersion in a sample of barred 
spiral galaxies. The bar strength is characterized by Qg, 
the maximal tangential perturbation associated with 
the bar, normalized by the mean axisymmetric force. 
It is derived from the galaxy potentials which are ob- 
tained using near-infrared images of the galaxies. How- 
ever, Qg is sensitive to bulge mass. Hence we also esti- 
mated bar strengths from the relative Fourier intensity 
amplitude {A2) of bars in near-infrared images. The 
central velocity dispersions were obtained from integral 
field spectroscopy observations of the velocity fields in 
the centers of these galaxies; it was normalized by the 
rotation curve amplitude obtained from HI line width 
for each galaxy. We found a correlation between bar 
strengths (both Qg and A2) and the normalized cen- 
tral velocity dispersions in our sample. This suggests 
that bars weaken as their central components become 
kinematically hotter. This may have important impli- 
cations for the secular evolution of barred galaxies. 
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1 INTRODUCTION 

In past decades there has been a large amount of dis- 
cussion regarding the secular evolution of bars in spiral 
galaxies. Numerical simulations show that gas inflow, 
star formation and the subsequent buildup of a central 
mass concentration (CMC) can considerably alter the 
shape of a bar and even result in the dissolution of the 
bar itself (Norman, Sellwood & Hasan 1996; Shen & 
Sellwood 2004; Kormendy & Kennicutt 2004; Hozumi 
& Hernquist 2005; Athanassoula, Lambert & Dehnen 
2005; Bournaud, Combes & Semelin 2005). However, 
the fraction of bars remains fairly high out to a redshift 
of Z'-^0.7 to 1 (Sheth et al. 2003; Elmegreen, Elmegreen 
& Hirst 2004; Jogee et al. 2004; Sheth et al. 2007) 
which suggests that either bars are long-lived or can re- 
form in galaxy disks. Simulations suggest that bars can 
dissolve but form again once the disk is cooled through 
gas accretion (Bournaud & Combes 2002) or through 
tidal interactions (Berentzen et al. 2004). 

Although there has been a considerable amount of 
numerical work on bar formation and dissolution, there 
is not much observational evidence. Molecular gas has 
a higher central concentration in barred galaxies com- 
pared to non-barred galaxies which suggests that bars 
are important for building up a CMC (Sakamoto et al. 
1999; Sheth et al. 2005). The bar length has been 
found to correlate with bar intensity contrast and cen- 
tral luminosity densities (Elmegreen et al. 2007). A 
more direct indicator of how bars change with mass 
concentration is the correlation of CMC and bar ellip- 
ticity (Das et al. 2003) which shows that bar ellipticity 
declines with increasing dynamical CMC. An enhanced 
CMC also results in the deepening of the central po- 
tential and gives rise to a dynamically hotter galaxy 
center. This will lead to a higher observed central stel- 
lar velocity dispersion in a galaxy (cr^). This effect has 
been seen in numerical simulations where a dynami- 
cally hotter center results when a bar dissolves (Fricdli 
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& Benz 1993; Hasan, Pfcnnigcr & Norman 1993) or 
weakens (Athanassoula. Lambert & Dehnen 2005). In 
this paper we investigate whether there is any obser- 
vational evidence for this heating effect by examining 
the correlation between the central velocity dispersion 
in bars and the bar strength. 

2 Central Velocity Dispersion 

The CT„ values used in this paper are from published 
integral-field spectroscopy (IFS) observations of the nu- 
clear regions of galaxies. We have thirty one galaxies 
in our sample (Table 1). The majority of ay values 
arc from the SAURON survey of nearby galaxies (de 
Zeeuw et al. 2002; Ganda ct al. 2006; Falcon-Barroso 
et al. 2006; Pclcticr et al. 2007); a smaller number 
are from the INTEGRAL and SPIRAL instruments on 
the William Hcrschel and Anglo-Autralian Telescopes 
respectively (Batchcldor et al. 2005) . The values for 
the remaining four galaxies are from observations using 
the GMOS instrument at the Gemini North Telescope 
(Barbosa et al. 2002). The velocity dispersion values 
that we have used are derived from two dimensional 
stellar velocity fields and not gas kinematics. 

One drawback of using these a^, values is that the 
aperture size is not uniform across the sample. Hence 
we have tried to bring the apertures to a common sys- 
tem using the aperture correction formula of Jorgensen, 
Franx & Kjaergaard (1995). The velocity dispersions 
were corrected to a radius Re/ 8, where Re is the ef- 
fective bulge radius of a galaxy. Thus (T„ was assumed 
equal to the average velocity dispersion within a radius 
Re/8 of the bulge of a galaxy. For five galaxies it was 
not possible to determine Re either because the bulge 
was not distinguishable or the image quality was poor. 
Also, since the galaxies span a wide variety of size and 
mass, we normalised the new aperture corrected veloc- 
ity dispersions (cTg) with the HI gas rotation velocity Vg 
for each galaxy. This also takes out some of the lumi- 
nosity/size effects for the galaxies. The Vg values were 
derived from the maximum HI rotation velocities (vh) 
using the galaxy axes ratios (q) to correct for inclina- 
tion (i.e. Vg = V}i / ^JX—q^ ); vh was obtained from 
the Hyperleda database. The ratio Oe/vg is thus an in- 
dicator of how kinematically hot the center is for each 
galaxy relative to its disk rotation speed. 

3 Bar Strength Derived from Near-IR Images 

There are several ways to quantify the strength of a 
bar in a galaxy. In this paper we have derived the 



bar strength in two ways; bar strength is assumed to 
be the maximum relative bar torque (Qg) derived from 
the gravitational potential of a galaxy. This is proba- 
bly the most robust estimate of bar strength but it is 
sensitive to the bulge mass or luminosity. The second 
way of estimating bar strength is to use the maximum 
of the relative intensity amplitude of the bar in the 
near-IR image (A2). Both methods require that the im- 
ages be deprojected before analysis. This assumes the 
disks are thin, but the presence of a less flattened bulge 
component could lead to artificial stretching of bulge 
isophotes. To minimize this effect, we decomposed the 
bulge from the disk, subtracted it from the image before 
deprojection, and added it back (assuming that bulge 
light is spherically symmetric) after dcprojccting the 
disk/bar light. The bulge components were separated 
from the disks using a two-dimensional multicomponcnt 
decomposition code which uses a Sersic model for the 
bulge, an exponential function for the disk, and either 
Ferrers' or Sersic's functions for the bar (Laurikainen, 
Salo, and Buta 2005). Effects of seeing were taken into 
account, using the values of the full width at half max- 
imum stored in image headers or provided in articles. 
The effective radius of the bulge (i?e) was estimated 
by integrating the fiux of the fitted bulge model. The 
images were obtained mainly from the 2MASS survey 
and some from previous studies. The fiters used were 
either K or H band. The bar parameters were derived 
as follows. 

(i) Qg : The gravitational potentials ($) were inferred 
from near-IR light distributions assuming that the light 
traces the mass. Bar induced tangential forces were 
calculated using a Polar method, as described in Lau- 
rikainen & Salo (2002) and Laurikainen, Salo & Buta 
(2004). In particular, the calculation applies an az- 
imuthal Fourier decomposition of intensity, including 
the even components up to m=20, which are then con- 
verted to the corresponding potential components (Salo 
et al. 1999). Two dimensional maps of the radial force 
{Fr) and tangential force {Ft) were calculated. The 
radial profile of the maximum tangential force at each 
distance is given by, 

^"^^^ ^ <l^.(.,0)l> 

where < |FR (r, </))| > denotes the azimuthally averaged 
axisymmetric force at each radius. The maximum in 
the Qt profile at the region of the bar then gives a single 
measure of bar strength, Qg. The main assumptions are 
that the mass-to-luminosity ratio is constant in the bar 
region, and that the vertical light distribution can be 
approximated by an exponential. The scale height, /i^. 
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was estimated from an empirical relation between hr/hz 
and the de Vaucouleurs type index T (de Grijs 1998), 
where hr is the radial scale length of the disk. The 
images were taken from the literature and are generally 
not very deep. Therefore, instead of estimating hr from 
the new decompositions, we used mainly the hr values 
from Baggett et al. (1998); if it was not available, we 
used other sources in the literature where deep optical 
images had been used to derive hr- For two galaxies hr 
was taken to be the mean value for our sample of 31 
galaxies. 

(ii) A2 : The same Fourier method gives us also the 
m~2 amplitudes of bar intensity contrast in the bar 
region. For some of the galaxies. Table 1 gives Qg and 
A2 but not the effective bulge radius. This is because 
application of the 2D decomposition method requires 
deeper images than the methods used to calculate Qg 
and A2. 

4 Estimating the Correlation 

Figure 1 shows bar strength Qg plotted against the nor- 
malized velocity dispersion o^jvg for 26 galaxies. Al- 
though there are 31 galaxies in the sample, we could 
use only 26 because the effective radius i?e could not 
be calculated for a few cases. The errors have been cal- 
culated using the standard error propagation equation 
and include the uncertainties in the observed quanti- 
ties. The majority of galaxies seem to follow a trend of 
decreasing bar strength with increasing central velocity 
dispersion. We quantified the correlation in two ways. 
The linear correlation coefficient for this sample of 26 
galaxies is r = —0.50 and the probabilty that they are 
from a random sample Pr is Pr < 1%. This method 
does not include the errors on both axes. A more ac- 
curate estimate would be a weighted correlation coeffi- 
cient, but this is difficult to obtain in practice (Feigelson 
& Babu 1992). Instead we used a simple Monte Carlo 
simulation that randomly samples the errors on both 
axes and determines a mean weighted correlation coef- 
ficient < r >. We used 50,000 linear fits and obtained 
a value of < r >= —0.46. The second method that we 
used to quantify the correlation was the Kendall- Tau 
coefficient which assigns relative ranks to the different 
values. This is perhaps a more robust way of examining 
the correlation especially when the sample size is rela- 
tively small as in our case. The Kendall- Tau coefficient 
for the 26 galaxies in Figure 1 is < tkt >= —0.35 and 
the probability that they are from a random sample is 
Pkt ~ 1.3%. 

Figure 2 shows the plot of A2 against the normalized 
velocity dispersion o-g for 25 galaxies. Here again we 
calculated the linear correlation coefficient which is r = 



—0.54 and Pr < 0.5%. When the errors are sampled 
using a Monte Carlo simulation we obtain a value of 
< r >= -0.51. The Kendall- Tau coefficient for the 25 
galaxies is < tkt >= —0.33 and Pkt < 2%. Thus 
both Figures 1 and 2 suggest that there is a correlation 
between the bar strength in galaxies and their central 
velocity dispersions. 

5 DISCUSSION 

The main result of this paper are shown in Figures 1 & 
2. There are 26 galaxies in the plots of which about half 
are intermediate type spirals and the remaining a mix- 
ture of early and late type spirals. Since the number 
of galaxies in each Hubble type is not very large, we 
cannot investigate trends within the different Hubble 
types. But from Figures 1 & 2, it appears that early 
type spirals have relatively lower central dispersions; 
this may be because they have larger bulges where ro- 
tational velocity is comparatively higher and the central 
velocity dispersion lower compared to the later Hubble 
types. 

Since the correlation is significant but not very 
strong we examined the two galaxies that define the 
higher and lower limits of Qg and A2. We looked at 
them closely to see if they are odd in some way, and not 
characteristic of the rest of the sample, (i) NGC 3162 
((Te/i'g=l-13) : This is an intermediate type spiral 
galaxy with a weak bar and prominent bulge. There 
may also be a ring in the center. Though the spiral 
arms are somewhat asymmetric, the nucleus is fairly 
undisturbed, (n) NGC 4314 {ae/vg=0.18) : This is a 
bright, early type galaxy with a strong bar and large 
bulge with a LINER type nucleus. There may be sig- 
nificant rotation in the nuclear region which may lower 
the (Te- Also, the mass may be more widely distributed 
over the bucle and hence not as centrally concentrated 
as in NGC 3162. There appears to be a ring of star 
formation in the nucleus as well (Gonzalez Delgado et 
al. 1997). Both galaxies are thus fairly normal and not 
unusually different from the rest of the sample galaxies. 

Figures 1 & 2 suggest that galaxies with dynamically 
hotter nuclei have weaker bars. It is now well estab- 
lished that a galaxy's central black hole mass and bulge 
velocity dispersion are correlated (Ferrarese & Merritt 
2000; Gebhardt et al. 2000). Later results show that 
the nuclear mass also correlates with the overall mass 
of a galaxy (Ferrarese et al. 2006; Hopkins et al. 2007). 
These results all indicate that the nuclear mass in a 
galaxy is intimately connected to the dynamics of its 
disk and halo. If so, then it is not suprising that ae in 
our sample of barred galaxies is correlated with Qg or 
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^2- It suggests that the growth of a central mass and 
evolution of the bar for spiral galaxies may be closely 
linked. 

Our results have important implications for the sec- 
ular evolution of barred galaxies. Simulations suggest 
that there may be several factors responsible for the 
dissolution of a bar. One is the CMC growth (Hasan & 
Norman 1990; Fricdli & Pfenniger 1991) which weakens 
the bar-supporting xi orbits and increases the fraction 
of chaotic orbits in the galaxy center. Second is the 
inflow of gas towards the galaxy center which results 
in the transfer of angular momentum to the bar wave 
which then weakens the bar itself (Bournaud, Combes 
& Semelin 2005). The buckling instability is an im- 
portant bar thickening mechanism that results in a 
boxy /peanut bulge which can temporarily weaken a bar 
(Raha et al. 1991; Berentzen et al. 1998; Athanas- 
soula & Misiriotis 2002; Martinez- Valpuesta, Shlosman 
& Heller 2006; Debattista et al. 2006). AU these ef- 
fects result in a more massive and dynamically hotter 
central component and a weaker bar. The correlations 
that we see in Figures 1 and 2 may be observational in- 
dications of this ongoing evolution. In particular, while 
the apparent drop of Qg with central velocity disper- 
sion might be an artifact caused by the bulge dilution 
effect (see Laurikainen, Salo & Buta 2004; a bias could 
follow since nuclear velocity dispersion and bulge mass 
are strongly correlated) , the similar correlation between 
the bar intensity contrast A2 suggest that the effect is 
real. The found correlation between A2 and <Je/vg is 
also consistent with Das et al. (2003) who found that 
the bar ellipticity (closely related to A2) drops with 
central mass concentration. 

The evolution of bars by secular processes in galaxies 
is an issue which is expected to gain more attention in 
the near future. Recent observational evidence shows 
that the fraction of strong bars in bright galaxies in- 
creases from under 10% at redshift z=0.84 to about 
30% in the local universe (Sheth et 2007). Also, it has 
been shown by (Laurikainen et al. 2007) that among 
the early-type barred galaxies the bulge-to-total flux 
ratios are on average smaller than in the non-barred 
galaxies. These results, together with ours, may indi- 
cate that bars evolve with their parent galaxies. 
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Normalized Central Velocity Dispersion (ct^) 

Fig. 1 The plot shows the bar strength Qg plotted against 
the normalized central velocity dispersion {(Je/vg) for 26 
galaxies from Table 1. The three galaxy types are plot- 
ted with different symbols; early- type spirals are SO/a to 
Sa (open circles), intermediate- type spirals are Sab to Sbc 
(filled triangles) and the late-type spirals are Sc to Sd (open 
triangles). The errors are marked on each axis. The dashed 
line is the best fit line and has the form y = —0.25a: -I- 0.39. 
Although there is some scatter, the overall trend suggests 
that bars weaken as their nuclei become dynamically hotter. 




Normalized Central Velocity Dispersion (aj) 

Fig. 2 The plot shows bar amplitude A2 measured from 
the near-IR images plotted against the normalized central 
velocity dispersion (a^/vg) for a sample of 25 galaxies. As 
in Figure 1 the open circles represent early type spirals, the 
filled triangles represent intermediate type spirals and the 
open triangles represent the late type spirals. The dashed 
line is the best fit line for the sample and has a form y=- 
0.45x + 0.75. 
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Table 1 Galaxy Data 
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0.391 


0.34 


0.62 


99 


2.4 




219.5 




(c) 


NGC 4293 


(R)SB(s)0/a 


early 


0.463 


0.36 


0.67 


42 


2.4 


4.65 


109.5 


0.35 


(c) 


NGC 4303 


SAB(rs)bc 


inter. 


0.861 


0.26 


0.44 


108 


2.7 


2.96 


66.3 


0.87 


(a) 


NGC 4314 


(Ri)SB(r'l)a 


early 


0.959 


0.44 


0.90 


43 


2.4 


5.46 


70.7 


0.18 


(c) 


NGC 4321 


SAB(s)bc 


inter. 


0.866 


0.18 


0.34 


101 


2.7 


8.62 


112.5 


0.45 


(a) 


NGC 4487 


SAB(rs)cd 


late 


0.659 


0.18 


0.18 


51 


2.4 


29.14 


86.7 


0.42 


(b) 


NGC 4593 


(R')SB(rs)ab 


inter. 


0.742 


0.31 


0.76 


105 


1.0 


4.32 


161.3 


0.44 


(d) 


NGC 4596 


SB(rs)0/a 


early 


0.716 


0.28 


0.67 


59 


2.4 


2.78 


92.5 


0.47 


(c) 


NGC 5005 


SAB(s)b 


inter. 


0.444 


0.15 


0.36 


203 


2.7 


3.03 


243.8 


0.79 


(a) 


NGC 5448 


(R;)SAB(rs)ab 


inter. 


0.454 


0.15 


0.65 


65 


2.4 


2.83 


188.4 


0.32 


(c) 


NGC 5585 


(R')SAB(s)d 


late 


0.646 


0.24 




42 


2.4 


0.17 


63.4 


0.59 


(b) 


NGC 5678 


SAB(rs)b 


inter. 


0.489 


0.18 


0.33 


103 


2.4 


2.43 


175.2 


0.54 


(b) 



"(i) deVA= de Vaucouleurs Atlas of Galaxies (Buta et al. 2007) 

''(ii) RC3= Third Reference Catalogue of Bright Galaxies (de Vaucouleurs, et al. 1991) 

^=(111) a=Batcheldor et al. (2005); b=Ganda et al. (2006); c=Peletier et al. (2007) d=Barbosa et al. (2006) 



